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The fragmentation functions of D°, D ± , D*, D*°, D* ± and A* at v^~10.6 GeV 
are measured with a data set of 102.7 fb^ 1 . Fragmentation model parametrizations 
(Peterson, Kartvelishvili, Collins-Spiller, Lund, and Bowler models) are compared to the 
data. The data at high x~l indicate a contribution of non-perturbative QCD processes. 
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1. Introduction 

The Belle experiment at the asymmetric e + e~ collided KEK-B in Tsukuba, 
Japan, is operating with an e~ beam energy of 8.0 GeV and an e + beam energy of 
3.5 GeV. The center-of-mass (cms) energy is adjusted to a y^^lO-58 GeV, corre- 
sponding to the mass of the T(45) resonance. At this cms energy, there is a resonant 
and a non-resonant component: the resonant T(45*) is produced with a cross sec- 
tion of <7~1.2 nb, the non-resonant continuum with <r~3.0 nb. The T(AS) decays 
into BB meson pairs with a branching fraction higher than 99%, and ~99% of all 
B{B) mesons decay into final states with hadrons containing charm quarks. In the 
continuum, the cross section of direct charm production is also high, i.e. er~1.3 nb. 
For about 10% of the time, data are recorded at y / s=10.52 GeV, containing only 
the non-resonant continuum component. 

Fragmentation is a QCD process which involves not only single partons, but a 
cascade of many partons, as it is shown schematically in Fig. 1. A single bare charm 
quark starts radiating gluons. The gluons themselves may a.) radiate additional 
gluons directly, or b.) convert into uu, dd or ss quark-antiquark pairs, while those 
quarks then might radiate additional gluons. If the distances between the quarks 
become too large (r>l fm), hadronisation starts and final state mesons and baryons 
are generated. 

The study of fragmentation has two aspects: On the one hand, as a theoretical 
aspect, QCD predicts hadronic confinement, but does not predict how an initial 
state quark as a pointlike object transforms into a final state hadron, i.e. an con- 
fined object with a formfactor. In most fragmentation models, the splitting of one 
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Parton Cascade 



Hadronizatloii 



Table 1. Decay modes used in the 
reconstruction. 




Df^TT+i^K+K-) 
A+~* p K-tt+ 

D*+-*D°Tr+(D"^K-n+) 

D*+^D+tt°Id+^K-tt+tt+) 

D*°^D it o (D ^K-it+) 



Fig. 1. Schematical description of the process e + e -^ce 
with subsequent fragmentation of the charm quark into a D* . 

parton to many partons (the "parton shower" ) is treated in perturbative QCD 
only, and hadronisation models are only phenomenological with minor predicting 
capabilities. Heavy quarks are in particular of interest, because their production is 
strongly suppressed in both the parton shower and the hadronisation. Thus, a final 
state charmed hadron contains a charm quark, which was produced with a high 
probability in the primary interaction. 

On the other hand, as a practical aspect, it is important to improve Monte-Carlo 
(MC) event generators. Several other analyses of direct charm production e.g. pro- 
duction of £>**, D s j etc. require MC simulations of charm fragmentation. In this 
paper PYTHIi^ 6.2 is used. The event generators use parameters, which have to 
be constrained by comparison to experimental data. 

The relevant variable for the investigation of the process dynamics is the momentum 
fraction x=phadron/ 'Pquark- As the quark momentum is a priori unknown and can 
only be approximated, in this paper x p =p/p ma x is used, where p m ax is corrected for 
the mass of the final state hadron, i.e. p ma x = \/s/4 — rn\ adron . The fragmentation 
function D{x p ) is then given as the normalized yield as a function of x p . Instead of 
the yield, related quantities such as the cross section da/dx p or BRxder/dSp can 
be used, where BR denotes the branching ratio into the particular final state. For 
the functional form of fragmentation functions, in PYTHIA the variable z is used 
instead of the variable x p , which is defined as z^(E+p z )hadron/(E+p z ) quar k. 
As will be shown below, in particular the high x region (i.e. x p ~l or z~l) might 
be sensitive to effects of non-perturbative G^CD. The x—1 limit is corresponding to 
the exclusive process e.g. e + e~^rcc— >D°D , in which on either side a bare charm 
quarks transforms into a hadron, which contains the identical charm quark as a 
valence quark, but additionally one light valence quark, light sea quarks and sea 
gluons. The interesting aspect is that, due to Phadron—Pquark, no momentum is 
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available to create the hadronic structure, but in fact the momentum balance is 
provided by the QCD vacuum. 

Fragmentation data of e + e~ — >cc at y/s~l0.6 GeV have been measured previously, 
but with lower statistical significance. More tha n 10 years ago, the ARGUS col- 
laboration and the CLEO collaboration publishecP^ data based upon integrated 
luminosities of 31.4 pb^ 1 and 35.8 respectively. Nowadays, equivalent data 

are recorded by BELLE in ~50 minutes. More recently, the CLEO collaboration 
published-^ D and D* fragm entatio n data based upon 8.9 fb^ 1 . D s and D* frag- 
mentation data were publishecPE21by CLEO and by BABAR based upon 4.7 fb^ 1 
and 23.4 respectively. 

2. Analysis 

A data set of 87.7 /6 _1 on-resonance and 15.0 off-resonance is used, and 

compared to Monte Carlo (MC) simulations corresponding to 217.0 fb" 1 . The MC 
simulations are based upon the event generator QQ9c0for e + e~— >cc, the Peterson^ 
fragmentation model, PYTHIA06.2 for the parton shower evolution and GEANlEl 
3.21 for the detector description. For specific simulations, also different fragmenta- 
tion models are used (see below). The decay modes used in the reconstruction are 
listed in Tab. 1. Charge conjugated modes are included. The efficiencies for particle 
identification of and are above 96% for all decay modes, the efficiency for 
protons is 81%. The particle misidentification probabilities are between 12% and 
26%. This paper represents a short version of the complete analysis which is pub- 
lished elsewhere^. With this unprecedented statistics it became possible to reduce 
the bin width in x p by a factor of ~3, and thus investigate the shape of the frag- 
mentation function with high precision, in particular in the high x p ~l region. In 
addition, fragmentation into the baryon A c is measured, for which a high statistics 
data set has not been available previously. 

Fig. 2 shows the mass distribution* 1 for the D° and the A+ (as an example for 
the meson and baryon reconstruction). The signals are shown for a momentum 
fraction 0.68<x p <0.70, which is close to the x p of peak position in the fragmentation 
function. The cross sections are calculated using the fitted signal yields. 
As mentioned above, at y / s=10.6 GeV there is a resonant and a non-resonant con- 
tribution. Thus, charmed hadrons can be produced by B decays (whereas the B 
mesons are produced by decay of the T(4S') resonance) or by fragmentation. In 
a first step, these two contributions have to be identified. Fig. 3 shows the cross 
section da/dxp as a function of x p . The down-left hatched histogram shows on- 
resonance data, the down-right hatched histogram shows off-resonance data. After 
normalisation using the integrated luminosities of the respective samples, the naive 
1/s dependence on the total hadronic cross section has been taken into account 
by a weighting factor (10.58 GeV/10.52 GeV) 2 . The contribution of B decays is 



a The mass distributions in Fig. 2 are s hown as mass difference distributions of the measured mass 
and the mass as given by the PDcffHI 
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dominant at lower a; p <0.5, as the B decay kinematics limits the maximum possi- 
ble momentum of the charmed hadron to about Phadron— 2-5 GeV/c, corresponding 
to x p ~0.5. The distributions in Fig. 3 are not yet efficiency corrected and only 
statistical errors are shown. 

3. Fragmentation Functions 

Fig. 4 shows the cross section da/dx p as a function of x p . For x p <0.5 only off- 
resonance data are used. For x p >0.5 the weighted average of off-resonance and 
on-resonance data is used. Fig. 4 is efficiency corrected, and represents the mea- 
sured fragmentation functions D(x), which can be compared to models. The inner 
error bars show the statistical, the outer error bars the total uncertainties. Fig. 4 is 
not yet corrected for feed-down from D* to D° . Details of the feed-down contribu- 
tion are described elsewhere^. The estimated systematic error from D** feeddown 
(e.g. £>£(2308) and D' x (2427)) is included and is <13%. Tab.Q]shows the functional 
form of the different fragm ent ation models whi ch a re considered. These are the 
Peterson^, the Kartvelishvil^l the Collins-Spillei^l the LuncP^ and the Bowlei^l 
models. The (1-x) term in all functional forms originates from the kinematics of 
radiation of one gluon. For every model a \ 2 is calculated from the bin-by-bin dif- 
ferences between the model and the data, while the model parameters are optimized 
in a fit procedure in order to find the optimum values. Fig. 5 shows the measured 
D* + fragmentation function in comparison to the models, and Tab. ^ lists the 
minimum / d.o.f.. The number of the degrees of freedom (d.o.f.) is given by the 
number of x p bins minus the free parameters of the fit (one or two, depending upon 
the model) . For the model comparison, the D* is chosen instead of the D° in order 
to avoid the feed-down correction. The x 2 values for all other hadrons can be found 
elsewhere^. The \ 2 is not based upon the pure functional form, but instead the 
MC generated fragmentation function according to the functional form. Thus, the 
detector resolution, as modeled by GEANT, as well as the parton shower, as mod- 
eled by PYTHIA, is taken into account. The models by Collins and Spiller and by 
Kartvelishvili are not included in the PYTHIA generator. Therefore, a reweighting 
procedure has been applied. For each event, the z and the p± values are stored, and 
the fragmentation function ansatz is recalculated. Tab. [21 shows the parameters at 
the minimum x 2 of the fit for all different hadrons. The fitted parameters indicate 
that the PYTHIA default parameters are not well suited for charm fragmentation. 
As an example, the PYTHIA default for the exponent a of the (l-x) a term in the 
Lund fragmentation model is a=0.30, but our fit gives as result a=0.68 at the min- 
imum x 2 - However, this is a global PYTHIA parameter, i.e. the parameter is valid 
for fragmentation of all quark flavors, and for e.g. the light quark sector the param- 
eter might be different from our result. As a result of the comparison, the Bowler 
and the Lund m odel s appear favored. This is consistent with an observation which 
was also reported^ for B meson fragmentation. The reason might by given by the 
exp(— brrij_/z) term with the explicit use of the transverse mass m± which appears 
in the functional form of these two models, but not in any of the other models. 
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4. High x region 

The measured data were used in a theoretical studjt^in order to investigate possible 
non-perturbative QCD effects. For this purpose, the fragmentation function is split 
into a perturbative and a non-perturbative part, i.e. D{x)=Dp(x)-\-Dmp{x), with 




using the normalisation N=J Q (1 — x) a x b . The factor /jvp indicates the relative 
fraction of the non-perturbative contribution and is determined by a fit. While for 
the cases of the D* mesons the non-perturbative contribution is /atp~21-25%, for 
the D° meson it turns out to be quite high with /np— 60%. The reason is, that 
a large fraction of the D° mesons arises from D* decay. The decay is very close 
to threshold, and thus the hadron momenta are almost equal, i.e. p(D*)~p(D a ). 
However, in the decay, a hadronic formfactor transition from a spin 1 hadron (i.e. the 
D*) to a spin hadron (i.e. the D°) occurs. The partons involved into this transition 
must be soft, and thus non-perturbative QCD power corrections are expected to 
occur. 

5. Summary 

Charm fragmentation functions were measured at GeV with high precision. 

The data were compared to five different fragmentation models, and the model 
parameters for the best fit to the data were given. The parameters can be used as 
input for PYTHIA. Among the five models, the Bowler and the Lund model are 
favoured. The data are posted^ to the DURHAM reaction database. 
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Fig. 2. Mass distributions for the D° (left) and the (right) for a momentum fraction 
0.68<x p <0.70. The histograms show the data, the dotted lines describe only the background, 
the full lines describe the signal and the background. 




Fig. 3. The cross section d<r/dx p as a function of x p for D° (left) and Ajt" (right). The down-left 
hatched histogram shows on-resonance data, the down-right hatched histogram off-resonance data. 




Fig. 4. Efficiency corrected cross section da/dx p as a function of x p . For rr p <0.5 only off-resonance 
data are used. For x p >0.5 the weighted average of off-resonance and on-resonance data is used. 

Table 1. Functional form of the fragmentation functions used in this analysis. The 
normalisation N is different for all functions. In the last column, the x 2 /d.o.f. from the 
comparison to the data for the case of the D* + (see Fig. 5) is shown as an example. 



Model 



Functional Form 



X 2 /d.o.f. for D*~* 



Bowler 
Lund 

Kartvelishvili 
Collins-Spiller N 
Peterson 



AT- 



- r ^(l- Z )' 2 exp| 

JV±(l-*)«exp(-^ 

Nz°"=(l - z) 
(^(l + ^l-I-^)- 



541.8 / 55 

965.6 / 55 
1271.1 / 54 
1540.7 / 54 
3003.0 / 54 
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Fig. 5. Efficiency corrected cross section BRxd<r/dx p as a function of x p in comparison to 
fragmentation models for D*+ fragmentation (black circles: MC, white circles: on-resonance data, 
white triangles: off-resonance data). Only statistical errors are shown. The parameters a, b of the 
models arc optimized by a fit procedure in order to find the best agreement between data and 
model. The deviations in the £~1 region are an indication for non-perturbative contributions. 



Table 2. The parameters of the fragmentation functions at the minimum of the x 2 /d.o.f. 
distributions, which can be used as input to PYTHIA. 







D° 


D+ 


Df 


A+ 


D*+ 


Bowler 


a\b 


0.12 | 0.74 


0.12 | 0.58 


0.12 | 0.68 


0.34 | 0.74 


0.22 | 0.56 


Lund 


a 


0.26 


0.45 


0.2 


0.55 


0.58 


Collins and Spiller 


4 


0.04 


0.055 


0.04 


0.04 


0.075 


Kartvelishvili 




4.6 


4 


5.6 


3.6 


5.6 


Peterson 


£c 


0.028 


0.039 


0.008 


0.011 


0.054 



